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ABSTRACT 

Kocevski, Ryde & Liang (2003, hereafter Paper I) proposed a semi-empirical 
fmiction (the KRL function) of gamma-ray burst (GRB) pulses, which could 
well describe those pulses comprising a fast rise and an exponential decay 
(FRED) phases. Meanwhile, a theoretical model which could give rise to this 
kind of pulse based on the Doppler effect of the expanding fireball surface was 
put froward in details in Qin (2002) and Qin et al. (2004, hereafter Paper II). 
To provide a physical explanation to the parameters of the KRL function, we 
try to fit light curves of the Doppler model (Qin 2002; Paper II) with the KRL 
function so that parameters in both models can be directly related. We pay 
our attention only to single GRB pulses whose profiles are that of FRED and 
hence employ the sample presented in Paper I (the KRL sample) to study 
this issue. We find from our analysis that, for light curves, which arise from 
exponential rise and exponential decay local pulses, of the Doppler model, 
the ratio of the rise index r to the decay index d, derived when fitted by the 
KRL function, increases quickly first, and then keeps nearly invariant with 
the relative width (relative to the timcscale of the initial fireball radius Rc/c) 
of local pulses when the width exceeds 2 (the relative width is dimensionless) . 
The rise and decay times of pulses arc found to be related with the Lorentz 
factor by a power law, where the power law index associated with the rise time 
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is less than that of the decay one and both of them are close to -2. In addition, 
the mean asymmetry shows a slightly trend of decreasing with lorentz factors. 
In plots of decay indexes versus asymmetry, there is a descending phase and 
after the phase there is a rising portion. We find that these long GRBs of 
the KRL sample are mainly associated with those light curves arising from 
CO- moving pulses with the relative width being larger than 0.1. Shown in 
our analysis, the effect of the co- moving pulse shape on the KRL hmction 
parameters of the resulting pulses is considerable and can be distinguished 
by the decay index d when the relative co-moving pulse width is less than 2 
(when the relative width is larger than 2, it would be difficult to discern the 
difference in the resulting pulse shapes). 

Key words: gamma-rays: bursts — methods: numerical analysis 

1 INTRODUCTION 

Light curves of gamma-ray bursts (GRBs) exhibit very complex and distinct morphologies, 
without any systematic temporal features (Fishman & Meegan 1995). A single spectrum, 
generally the Band function spectrum, seems to be a universal characteristic of GRBs (Piran 
et al. 1997). The highly variable temporal structure of GRBs may provide the clue to disclose 
the puzzle about the spectrum observed (Sari et al. 1997). It was previously found that some 
individual peaks decay more gradually than they rise (e.g. Fishman et al., 1994; Fenimore 
E. E., 1999). Many GRBs can well be decomposed into fundamental pulses (Norris et al. 
1996; Lee, Bloom, & Petrosian 2000a, 2000b), which comprise a fast rise and an exponential 
decay (FRED) phases in general. The pulse-shape light curves as the elementary events of 
GRBs probably provide the intrinsical information. Time asymmetry in GRB light curves 
has been discussed by several authors (see Barat ct al. 1984; Norris et al. 1986; Kouveliotou 
ct al. 1992; Link, Epstein, & Priedhorsky 1993; Nemiroff et al. 1994). The rise time (tr), 
decay time {td) , and full width at half maximum {FWHM) are the main factors concerned 
(McBreen et al. 2002; Ryde et al. 2003; Kocevski, Ryde & Liang 2003, hereafter Paper I). 
Norris et al. (1996) utihzed a single flexible empirical function to represent pulses in long 
bright GRBs, and found the ratio of rise-to-decay times of pulses to be unity or less. Shown 
in Ryde & Svcnsson (2000, 2002) are some other fitting profiles. 

To fully model the FRED light curves or individual pulse shapes, Kocevski et al. (2003) 
(Paper I) put forward an empirical function (the KRL function) described by 
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where tm is the time of the maximum energy flux, Fm, of the pulse. The shapes of single 
GRB pulses are confined by the two parameters r and d. It is a very functional and flexible 
model based on the physical first principle and the well-established empirical correlation 
between Ep and flux. This function could be adopted to characterize some GRB pulses such 
as FRED hght curves. However, what can parameters r and d reveal if the FRED pulses 
described by the KRL function arise from fireball sources? 

Recently, a theoretical pulse model (called the Doppler model) was proposed (see sec- 
tion 2 below) to relate the observed characteristics of GRB pulses with parameters of the 
intrinsic radiation in the framework of fireballs, where the Doppler effect associated with 
the expanding fireball surface is the key factor to be concerned (Qin 2002; Qin et al. 2004, 
hereafter Paper 11). Illustrated in this model, the profile of an observed pulse would mainly 
be determined by the intensity of co-moving pulses, I{Tg) (see Paper II). It was revealed that 
most of the resulting pulses possess a characteristic of FRED, even though the co-moving 
pulses concerned are diverse significantly. 

Hinted by the common FRED feature shown in the two models, a primary goal of this 
paper is born. Wc want to know how parameters r and d would be associated with intrinsic 
quantities such as the Lorentz factor P and the co- moving pulse width Arg (sec below), once 
co-moving pulses are provided (see Paper II). Some prevenient investigations did't take into 
account the infiuence of co-moving pulses on the observed light curve (see, e.g., Ryde 2004). 
However, as mentioned above, it was shown in the Doppler model that the observed pulse 
shape (especially its peakedness) is mainly associated with the co-moving pulse shape and 
width. 

In the following, we will investigate the impacts of exponential rise and exponential decay 
co-moving pulses on the observed light curve in detail. In section 2 we will introduce formulas 
of the Doppler model to describe the observed FRED pulses. In section 3, we will give a 
reason for choosing the types of co-moving pulse and will simulate a lot of resulting pulses 
and then will fit them with equation (1), and in doing so, some parameters associated with 
the observed light curve will be presented. In section 4, we will contrast our results with 
those derived from observation. Our results will be discussed and summarized in section 5. 
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2 FORMULAS OF THE DOPPLER MODEL EMPLOYED IN THIS PAPER 

We list in this section basic formulas, which were derived and presented in Paper II, of 
the Doppler model so that it would be convenient for us to employ or refer to them in the 
following analysis. 

Adopting the same symbols assigned and used in Paper II, the count rate equation 
derived from the expected flux emitted from an expanding fireball surface is expressed as 
(see Paper II) 

J I{re){l + PTef{l-T + Te)dTer-'^dv 

^^^^ ^ ^'-^ r3(l-mi + A-)^ ' ' 

where (5 — y/V^ — l/F, F is the Lorentz factor, Cq is a constant which includes the luminosity 
distance D between the fireball and the observer and other factors (see Paper II), /(r^) 
represents the development of the intensity of the co-moving emission (in this paper, the co- 
moving pulse), go^u{,T-'o,0) describes the rest-frame radiation mechanism, and 1^0,61 is the rest 
frame emission frequency which is related with the observed frequency v by the Doppler 
effect. In equation (2), the count rate concerned is defined within an energy channel of [i^i, 
i/2]. As shown in Paper II, letter "i" denotes an absolute value of time which is normally 
defined (it is always in units of s), and the Greek letter "r" represents a relative time scale 
(in units of 1) corresponding to "f . The two types of quantity are related by (see Paper II) 
_ t - tc - D/c + Rc/c 



Rc/c 



(3) 



Te,min = , 

and 

t0,max tc /r:\ 

'^9, max = ^ 1 ) [P) 

ricjC 

where tc and Rc are constants, and we assign te^mm < ^6> < td,max (^-nd hence Tg,rnm < 
Te < TQ^rnax)- The limits of the integral of tq in formula (2) are fe^rnin = ?Tiax{re^mj„, (r — 

1 + COs6'maa;)/(l - (3 COS 6 „iax)} and Te,max = min{T0^rnax, {t - 1 + COs9min){l - P COS 9min)} , 

where Omin and Omax are the lower and upper limits of angle 6, respectively (see Paper II). 
One could find from these definitions that r and tq are dimensionless quantities. 
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Note that, the value of could be arbitrarily chosen (it depends on the time we assign 
to that moment). For the sake of simplicity, we take 

t, = {R, - D)/c (7) 

in the following analysis, where Rc is the initial radius of the fireball (measured at time tc). 
In this way, one gets 

Substituting this relation into formula (2), the count rate as a function of observation time 
t could then be plainly illustrated, so long as Rc is available (or assumed). 

In this paper, we take the generally adopted spectral form proposed by Band et al. 
(1993), the so-called Band function 



{f-tr-^'>expWo - «o)(^)^+* > f^) ' 



2+ao ' /^V/^u ^'^VQ,p' ^'^o,p 2+ao 

as the rest frame spectrum g'o,i/(^o,6»)- 

One might observe that the profile of the KRL function are determined by parameters 
r and d. However, in the Doppler model, the characteristics of an observed pulse would 
be obviously influenced by the Lorentz factor and the co-moving intensity of radiation (see 
Paper II). Thus, by fitting light curves of formula (2) with function (1), one might be able 
to tell how parameters r and d and various widths of the KRL function are related with the 
Lorentz factor and co-moving pulses, assuming that the sources observed are undergoing a 
fireball stage. 



3 GENERAL ANALYSIS ON THE PULSES OF FIREBALL SOURCES 

In this section, we study the light curve of formula (2) in a particular case: the adopted co- 
moving pulse is an exponential rise and an exponential decay one. Various sets of parameters 
r and d of the KRL function would be obtained when fitting light curves of formula (2) 
associated with different Lorentz factors and different widths of the co-moving pulse with 
equation (1). In this process, we consider the radiation from the whole spherical surface of 
the fireball, although the contribution from the area of ^ > l/F is insignificant. (Impacts of 
different emission areas in the fireball surface on the resulting pulses were shown in Paper 
II.) Thus, we adopt = and 9.„iax = 7r/2. Not losing the generality, we assign Cq = 1. 
In addition, we take i/i — BOi/Qp and 1/2 — lOOi/Qp, where is the break frequency of the 
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rest frame Band function spectrum, so that when assigning i/o,p = lkeVh~^ the energy range 
would correspond to the second channel of BATSE. 

It is supposed that a very violent coUision of two shells gives rise to the observed GRB 
pulse behavior, and the increase (or decrease) of co-moving emission is proportional to the 
total radiation in the rise (or decay) phase of the co-moving pulse, namely 



where Jq and Tg^ are the integral constants. Associated with the well-known mechanisms, 
the rising part of this co-moving pulse might connect to the shell crossing time and the 
decay portion might relate to the cooling time. (It is interesting that a pulse form with an 
exponential rise and an exponential decay phases was previously adopted as an empirical 
function to describe some observed GRB pulses by Norris et al. 1996). The spiky form of 
this co-moving pulse implies that the physical interaction processes of shells are greatly 
impetuous and rapid. 

In the process of the shell collision, two intrinsic timescales, the cooling timescale of 
electrons and the shell crossing timescale, should be considered simultaneously, and they 
together would cause the width of the co-moving pulse. In terms of the Doppler model, 
co-moving pulses would be significantly altered by the expanding fireball surface and then 
would lead to the observed forms of pulses (see Paper 11). An alternative interpretation to 
this is that the two timescales together with the curvature timescale which is due to the 
relativistic kinematics of expanding shell would give rise to the observed pulses. In many 
cases, the curvature time and the shell crossing time dominate over the cooling time (See, 
Ryde et al. 2003). Spada et al. (2000) suggested that this case will occur at a distance 
i? < 5 X lO^^cm, but for larger radii the radiative cooling time will become the dominant 
contribution to the pulse duration (in the first case the cooling timescale would be relatively 
small while in the last case it would be relatively large). 

For convenience, the width of co- moving pulse (11) is defined as 

— 'Te,max — Te,min- (12) 

Note that the quantity Arg determined by eqs. (5) and (6) is also dimensionless. In the 



I —-^I{Tg)dTg {ill the decay phase) 
where cr^ and cr^ are two positive constants. Integrating this equation leads to 




(10) 




(11) 
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following analysis we assign Tg^min = ^o'n Te,max = o'd and tq^q = in order to check whether 
the ratio of the rise time to the decay time of co-moving pulses can influence the observed 
light curves. Here we consider two kinds of spiky co-moving pulses. One concerns (7^ = 
(case 1) and the other is associated with (7^ = 4(7d (case 2). 

3.1 In the case of cr^ = aa 

Here, we consider co-moving pulse (11) with cTj. = ad, which implies that the cooling time 
and the shell crossing time are comparable. 

3.1.1 Impact of the co-moving pulse width 

To investigate how parameters r and d are related with the co-moving pulse width, we 
calculate various light curves of formula (2), arising from co-moving pulse (11) with Arg 
= 0.001, 0.002, 0.005, 0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1, 2, 5, 10 respectively (The reason to 
choose these values is discussed in Appendix B.), and then fit them with the KRL function. 
For the rest frame Band function spectrum we take ao — —1 and (3o — —2.25, and for 
the Lorentz factor we consider F = 2, 5, 10, 100, 1000, and 10000. Parameters r and d are 
directly obtained from the fit. The rising width tr and the decaying width td of the observed 
pulse would be obtained from equation (1) when the fitting parameters are applied. (Note 
that, quantities tr and td are consistent with those denoted respectively in Paper I.) We 
define the asymmetry of the observed pulse as tr/td, the ratio of the rise fraction timescale 
of the pulse to the decay fraction timescale. 

(i) Relations between these parameters (r, d, r/d, tr, td, tr/td) and the co-moving pulse 
width At0 are shown in Fig. 1. Shown in panel (a), there lie minimums of d between Arg — 
0.03 and 0.1. From panels (a), (b), (c) and (f), we find that d, r, r/d and tr/td are independent 
of Arg when Arg is larger than 2. In terms of the KRL function, the result suggests that 
profiles of the observed light curves would not be distinguishable when the width of co- 
moving pulses is large enough (say Arg > 2), which is in agreement with what revealed in 
Paper H (see Paper H Fig. 3). According to equation (4), a large value of Arg suggests a 
relatively small value of the radius of the fireball, and this in turn indicates that the curvature 
timescale is relatively small. In this situation, it would be reasonable when the cooling time 
plus the shell- crossing time dominate far over the curvature timescale. However, as shown in 
panel (f), the pulse asymmetry would decease rapidly with the decreasing of the co-moving 
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Figure 1. Relations of parameters (d, r, tr, td, r/d, and tr/td) of the observed light curves and the co-moving pulse width 
(Arg) for different lorentz factors (F = 2, 5, 10, 100, 1000, 10000) respectively. Symbols corresponding to different values of T 
are displayed in panel (a). 

pulse width when Ar^ < 2. Panels (d) and (e) indicate that increases with the increasing 
of Are at all time, while td is sensitive to Ar^ only when Ar^ is large enough (say Arg > 0.1). 
This must be due to the fact that when the curvature timescale dominates far over the two 
other timescales, the decay phase of the light curve would be determined by the curvature 
effect and its profile would remain unchanged when Ar^ < 0.1 (approaching the so-called 
standard form defined in Paper II; see also Paper II Fig. 3). A contrast between panels (c) 
and (f) suggests that the ratio of r to d doesn't adapt to characterize the asymmetry of the 
observed pulses. 

(ii) Displayed in Fig. 2 are the relations between parameters rf, r, r/rf, tr/^cz? FWHM, 
and tr/FWHM. Shown in panel (a), the smallest value of index d could be detected in the 
range of 0.5 < r < 1. Panels (b) and (f) also demonstrate the existence of a minimum of 
d, with the minimum corresponding to smaller values of tr/t^ and t^/ FWHM when the 
Lorentz factor becomes larger. One finds from panel (e) that, for a given Lorentz factor, 
the pulse asymmetry is very sensitive to the co-moving pulse width when the latter is very 
small, and when the latter becomes large enough, the asymmetry would become invariant. 
This is in good agreement with what illustrated above. In addition, panel (e) shows that the 



Physical implication of the KRL pulse function of gamma-ray bursts^ 9 




FWHM FWHM U 



Figure 2. Relations between parameters d, r, r/d, tr/tci, FWHM, and U/FWHM, where we adopt T = 2, 5, 10, 100, 1000, 
10000 respectively. Symbols associated with different Lorentz factors are denoted in panel (a). 

larger the Lorentz factor the narrower the pulse observed, which is in consistent with what 
previously known (see, e.g., Fenimore et al. 1993). 

We find that characteristics of the relationships displayed in Figs. 1 and 2 are the same 
for different Lorentz factors. 

3.1.2 Impact of the Lorentz factor 

Here, we study how parameters d, r, r/d, tr, td and tr/td are related with F when different 
values of Ar^ are adopted. We take F = 2, 5, 10, 70, 100, 150, 200, 500, 1000, 10000 
respectively. Five typical values of the co-moving pulse widths, Ar^ = 0.001, 0.01, 0.1, 1, 
and 10, are adopted. The indexes of the rest frame Band function spectrum are the same 
as those adopted above. In the same way, light curves of (2) associated with various sets of 
these intrinsic parameters are calculated and then are fitted with equation (1). 

Shown in Fig. 3 are the relations of parameters of the KRL function associated with 
the observed light curves and the Lorentz factor. As shown in panel (f), the effects of F on 
the pulse asymmetry are negligible, although a weak anti-correlation between t^/^d and F is 
visible. The following relations could be concluded from panels (d) and (e): 
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Figure 3. Relations of parameters (d, r, tr, ^d, r/c/, tr/^d) observed light curves and the Lorentz factor F associated 

with different values of the co-moving pulse width (Arg = 0.001, 0.01, 0.1, 1, 10) respectively. Symbols are denoted in panel 
(a). 

tr OC (13) 

and 

U OC T'^^. (14) 

Values of and ad could be obtained by performing a linear fit to the logarithmic format 
data of the two panels. The results are listed in Table 1. The relation Atfwhm oc 
was obtained by Qin et al. (2004) (Paper II) when an extremely narrow co-moving pulse is 
concerned. The index of —2 is nothing but merely a result of time compression effect caused 
by the forward motion of the ejecta (see Appendix A). The values of and presented in 
Table 1 are close to —2, which must be due to the same time compression effect. Note that 
the observed timescale lasts a little longer (as the indexes are slightly larger than —2) than 
what the time compression effect suggests. This could be understood when one recalls that 
what we consider in this paper is a fireball surface rather than an ejecta moving towards the 
observer (in the latter case, ^ = 0) and the emission of the surface lasts an interval of time 
(if the emission is extremely short, we come to Paper II equation [44]). Shown in Table 1 we 
observe that ad is slightly larger than a^. This suggests that the rising part of the observed 
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Table 1. Power law indices Or and o^j derived for cases 1 and 2 
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0.001 


-1.95 


2.7 


-1.92 


3.1 




0.01 


-1.98 


3.0 


-1.93 


2.6 


0"r = 0"(i 


0.1 
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-1.99 
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2.8 
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-1.99 
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Figure 4. Relations between parameters of the observed light curves associated with Arg = 0.001, 0.01, 0.1, 1, 10 respectively. 
Symbol are denoted in panel (a). 

pulse is less affected by the fireball surface than the decay portion is (the less affected by 
the fireball surface the closer the index to —2). 

(ii) Relations between parameters shown in Fig. 2, the characters of the observed light 
curves (ie., d, r, r/d, tr/td, FWHM, and tr/FWHM), are displayed in Fig. 4, where different 
data points correspond to different values of F (data points associated with the same value 
of Arg are denoted by the same symbol). 

Panel (e) shows that, for very narrow co-moving pulses, asymmetry of the observed 
pulses would keep to be invariant with FWHM, while for wider co-moving pulses, the 
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asymmetry would decrease first and then would keep to be invariant with the decreasing of 
FWHM (here, since the co-moving pulse width is fixed, the decreasing of FWHM would 
be caused by the increasing of the Lorentz factor). Panels (b), (c) and (f) are similar to 
the corresponding panels in Fig. 2. They imphes that: (1) there is a correlation between 
parameters d and tr/FWHM in tr/FWHM > 0.2 (or Atq > 0.1) and an anti-correlation 
between the two quantities in tr/FWHM < 0.2 (or Atq < 0.1) (see section 4 for a detailed 
discussion), no matter what the Lorentz factors are. (2) a power law relation probably exists 
between parameters r and tr/td- Let us assume 

tr/td oc r^ (15) 

Fitting the logarithmic format data with a linear function, we obtain e = 0.578. It shows 
that index r and the pulse asymmetry tie in, and the former can measure the latter as long 
as the co-moving pulse is known. 

The result that the pulse asymmetry for long bursts will decrease as the FWHM narrows 
supports previous conclusions found by Norris et al. (1996) (see also e.g. Reichart et al. 2001). 
As shown in Paper II, a small value of FWHM could be caused by a small co-moving pulse 
width, or a small fireball radius, or a large Lorentz factor. Of the three factors, the third is 
the most sensitive one according to the Doppler model. 



3.2 In the case of cr^ = 4(7^ 

Now we study the case of adopting co- moving pulse (11) with (7^ = 4(7^ (case 2), which 
corresponds to a relatively fast coohng timescale of electrons. 

In order to find out whether there is any difference between the two cases, we perform 
the same analysis as that in case 1. The only difference in this analysis is that we replace 
ar = Od with Or = 4:ad for co-moving pulse (11). It is surprised that characteristics shown 
in the corresponding relations in the two cases are very similar. Conclusions drawn from 
case 1 holds in case 2. However, parameter d shows a difference in the two cases, which is 
illustrated in Fig. 5 (presented in this figure is also the comparison between the values of r 
for the two cases). To produce the data of this figure, we take F = 100. As suggested by the 
figure, index d is more sensitive to the co-moving pulse than r is. This enables us to relate 
an observed pulse with the corresponding co-moving pulse by parameter d so long as Arg is 
smaller than 2. It suggests that a lager value of d may correspond to a faster cooling process. 
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Figure 5. Relations between parameters (d and r) of the observed light curve and the co-moving pulse width for cases 1 and 2. 
Here we take F = 100. Relations associated with d and r are presented in panels (a) and (b), respectively. Symbols are denoted 
in panel (a). 

Once Arg becomes larger than 2, it would be difficult to discern different co-moving pulses 
from d. 

The analysis shows that parameters ar and ad in case 2 are larger than those in case 1. 
This indicates that the influence of the co-moving pulse in case 2 on observations is greater 
than that in case 1. 

In case 2, regarding relation (15) we get e = 0.464 by a linear fit. It shows that the same 
conclusion obtained in case 1 holds in this case. 

4 COMPARISON WITH THE OBSERVED DATA 

Let us contrast Fig. 4 panel (f) with Paper I Fig. 12 to provide a direct comparison between 
the observed data and the expectation of the Doppler model in the relationship between d 
and tr/FWHM. In Fig. 4 panel (f), data in the decaying portion of the relationship curve 
(the data associated with Atq < 0.1) are defined as sample 1, while those in the rising 
portion (the data associated with Atq > 0.1) are called sample 2. The corresponding sets of 
data in case 2 are also called sample 1 (the data associated with Arg < 0.1) and sample 2 
(the data associated with Arg > 0.1) respectively. The data presented in Paper I include 77 
individual pulses with time profiles longer than 2 seconds, which we define as sample 3. 

One might observe that analysis performed in the previous sections is based on the 
concept of r which is dimensionless. Relation between this quantity and the observed time 
t is shown in equation (8). As r is proportional to t, parameter tr/FWHM would be the 
same in both definitions. Fortunately, d is dimensionless. We therefore can directly compare 
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Figure 6. Comparison between sample 3 and samples 1 and 2 in the planes of d vs. tr/FWHM and number vs. d, where 1) 
panel (a) is associated with case 1 while panel (b) corresponds to case 2. [Note that, data of samples 1 and 2 plotted in panel 
(a) are merely the replicate of Fig. 4 panel (f).] The vertical line denotes the position of tr/FWHM ~ 0.2 (or Arg ~ 0.1) which 
divides data of sample 1 from those of sample 2. Symbols are denoted in panel (a). 2) Panels (c) and (d) shows distributions 
of index d for case 1 and case 2 respectively. The assumed lognornmal distribution of Arg [see eq.(16)] yields the theoretical d 
distribution as the lorentz factor is taken one typical value, say, P = 100. Same symbol meanings as in panel(c). 

data of samples 1 and 2 with those of sample 3. Plots of d vs. tr/FWHM for cases 1 and 
2 together with the data of sample 3 are shown in Fig. 6 panels (a) and (b) respectively. 
Presented in the figure is also the division (at tr/FWHM ~ 0.2) between samples 1 and 2. 

We find in surprise that sample 3 is well within the range of sample 2, while it is com- 
pletely irrelevant to sample 1. In addition, sample 3 shows a positive correlation between d 
and tr/FWHM just as what sample 2 shows. To investigate whether samples 2 and 3 have 
indeed the same distribution on a certain significance level, we take a general K-S test to 
the two-dimensional distributions of the two samples. We adopt the effective D definition 
of two-dimensional K-S statistic D as the average of two values obtained by above samples 
individually, (see. Press et al.(1992)) The K-S statistic and the indicated probability (signifi- 
cance level) are listed in Table 2. Above results in both case 1 and case 2 imply that sample 
2 and sample 3 are not significantly different in terms of statistics. Based on the fact that 
sample 3 is composed of 77 individual pulses with durations longer than 2 seconds, we thus 
deduce that long GRBs could be mainly superposed by wider rest-frame radiation pulses 
with At0 > 0.1. Motivated by the viewpoint that the short bursts with Tgo < 2.6s have 
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Table 2. Parameters gained by two- and one-dimensional K-S test for respective cases 1 and 2 



Panels 


KS 


Probability 


Degrees of freedom 


a 


0.3075 


0.0475 


nl=30, n2=76 


b 


0.4145 


0.0022 


nl=30, n2=76 


c 


0.2237 


0.0376 


nl=n2=76 


d 


0.2895 


0.0026 


nl=n2=76 



different temporal behaviors compared with the long ones and they may actually constitute 
a different class of GRBs (e.g. Norris, Scargle, & Bonnell 2001), we suppose sample 1 with 
Are < 0.1 might hold the characteristic of short GRBs. The fact that sample 2 associated 
with At0 > 0.1 is consistent with sample 3 in distributions seems to demonstrate at least 
some of sources in sample 3 could be described by Qin's theoretical model in a sense of 
practice. 

To check if the distribution of index d in sample 3 is indeed expectable, we plot the 
distribution of d that is found from Fig. la for a assumed distribution of Ar^. The rise times 
(t^), fall times (t^), FWHM, as well as the time intervals between pulses had been measured 
and found to be consistent with lognormal distributions for both short and long GRBs 
(jMcBreen et al. 2002|) . Now, let us also consider Gaussian distribution of the logarithmic 
format of quantity Ar^: 

pilog^r,) « exp|-i(^^«<^I^±H5)^|, (16) 

So, many random values of Ar^ could be yielded so long as the reasonable range of Atq 
can well be determined. In Fig. la, suppose we take F = 100, a correlation between d and 
Are within the range of 0.1 < Arg < 2 (Note that, resulting pulses within this range are 
correlated with sample 3 and can then be distinguished with parameter d in KRL function) 
offers us an opportunity to gain the index d corresponding with Ar^. Plots of distribution 
of d for cases 1 and 2 together with the data of sample 3 are shown in Fig. 6 panels (c) 
and (d) respectively. Likewise, the K-S test to above distributions is also made to give the 
statistic and indicated probability listed in Table 2. We find the two distributions are not 
significantly different from viewpoint of statistics. 



5 CONCLUSIONS AND DISCUSSIONS 

According to above-mentioned analyses, we can draw the underlying conclusions: 

First, effects of different sorts of spikier co- moving pulses on the observed pulse shapes 



16 Zhang and Qin 

could be distinguished by shape parameter d provided Atq is less than 2 or so. However, 
parameter r in KRL function and pulse asymmetry could be bound with the relationship 
tr/td oc r^ Secondly, the asymmetry of observed pulses exhibits a slight trend of decrease 
with the increasing of F, whereas it increases quickly with the increasing of Atq when 
< 2, beyond this range, it will keep invariant with Atq. In observer framework, we find 
the asymmetry increase quickly first with FWHM and then behaves nearly independent of 
FWHM for an assumed F or a lager Atq. Thirdly, two power law relations, ~ F"'' and 
td ~ F"'* have been surprisingly found to show the width of observed pulse is highly sensitive 
to lorentz factors. Along with previous result (Paper II eq.[4]), we attribute these properties 
to the so-called time compression effect, which is purely kinematical and independent of 
physical mechanism, as long as the emission comes from an expanding fireball. Further, the 
difference of indexes and in case 1 from those in case 2 suggests that diverse intrinsic 
emission processes may cause distinct influences on observed profiles. 

Following from the discussion in §4, for reasons not presently understood, the result 
that there is a decay portion (sample 1) in plot of d vs. tr/ FWHM which disappears in 
Kocevski's plot (Paper I, Fig. 12) is rather surprising. We believe that there might be some 
bursts arising from short co-moving pulses associated with Atq < 0.1 and their corresponding 
pulses would be observed if our sample is large enough, and in this case the corresponding 
data would be located within the descending portion in this plot. The two sub-classes of 
GRBs show completely different behaviors though they can be explained with the same 
emission mechanism, which in general being ascribed to synchro-Compton radiation via 
internal shocks, not the external shocks (see, Piran et al. 1997, Nakar et al. 2001 & McBreen 
et al. 2001, 2002). Based on this comprehension, we infer that the decay portion probably 
encompasses the short bursts or at least some of them. 

As mentioned above, our whole investigations are built on the assumptions of the fireball- 
model and the isotropic radiation. Further, all the co-moving pulses concerned in this paper 
have a certain width, that is Arg > 0. Previous studies on light curves of a co-moving S- 
function pulse (see. Paper II) found that this very narrow pulse will lead to a standard form 
light curve, which has only decay phase in the resulting pulses due to pure curvature effect. 
As it shows in Fig. 1 that the difference in contribution of Atq to rise times and decay times 
of light curves is visible. This shows the rise phase of pulses as a result of the contribution 
of Atq would reflect not only the energizing of the shell but also the radiative cooling of 
electrons, while the decay portion of the observed pulse could be fully characterized by 
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all above-mentioned timescales (namely, curvature time , shell-crossing time and radiative 
cooling time). In the case of Arg ^0.1, the Doppler model is not the major contributor of 
the pulse shape and indeed the co- moving behavior could be important. The conclusion is 
just in excellent agreement with that of Spada et al.(2000) 

By simulating many resulting pulses with several sorts of co-moving pulses, We find 
the resulting pulses modelled by Qin's model are mainly determined by their corresponding 
co-moving forms and lorentz factors. Assuming one monotonic rise (or decay) function is 
selected to stand for the co-moving form, the resulting pulses will show concave (or convex) 
in phase of rise. In addition, we find the kurtosis especially the rise part of resulting pulses 
mainly originates from the spiky co-moving forms as long as Atq is not very narrow enough, 
otherwise, these peaked resulting pulses will behave greatly similar to the so-called standard 
forms. In other cases, we can always achieve the flat-peaked resulting pulses provided that 
the decay phase of co- moving pulses is in existence. 

From an viewpoint of observation, how to choose the co-moving pulse form for an ob- 
served GRB is an urgent issue. Within the range of 0.001 < Are < 2 (namely, 0.5 s <"t" 
<1000 s), the shape of resulting pulses could be well distinguished by shape-related pa- 
rameter d. This may enables us to regard d as a probe to speculate on what the detailed 
co-moving forms are. On the other hand, we also need synthetically take into account the 
special physical process. Motivated by these considerations, we might give the relatively 
correct co-moving pulse form which is then utilized to fit to observed data of GRBs. With 
the best fit-of-goodness, the likely parameters in rest-frame such as F, Atq and Rc could be 
derived. Once the co-moving pulse shape is adopted, it can be applied not only to derive 
some parameters in co-moving framework, but also to constrain physical emission mech- 
anism and improving other theoretical models. In this paper, we primarily focus on the 
theoretic analysis of pulses. Some results could be primary and need to be approved by 
more observations. 
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APPENDIX A: TIME COMPRESSION EFFECT 

Here we show how the index of the power law relation between an observed timescale and 
the Lorentz factor is —2. 

Let us consider an ejecta moving towards the observer with a velocity of v = cjS, which 
emits two photons at different times. Suppose the ejecta emits the first photon from distant 
D at its co-moving time tg^i, and emits the second one from distant D — v{tg^2 — te,!} at 
time tQ2- It is obvious that the observer receives the two photons at its observed time 
ti = te,i + D jc and t2 = te,2 + [D — v{tg^2 — ^6»,i)]/c, respectively. This leads to 

t2-ti= t0,2 - te,i - (3{te,2 - te,i). (Al) 
Applying 1 — /3 ~ l/2r^ one gets 

t2 - ti - ^j,, . (A2) 

APPENDIX B: DISCUSSION ABOUT TAKING THE VALUE OF 
CO-MOVING PULSE WIDTH(Are) 

Since the fireball becomes optically thin when radiation of 7-ray begins, the size of it can 
be estimated, in general, to be 10^^^^ ''cm (Piran, 1999). Thus 

Rjc ~ lO^'^sec (Bl) 
From (4), we can get 

By using another relation 

At, = ^ f , (B3) 
1 — p cos 



(Qin et al. 2004) 



we can deduce 

At 

= Rjc{l-f3cos9) ^^'^ 
Norris et al.(1993, 1996) had pointed the entire range of all pulse widths is ~ 10ms — 2s, 
that is to say 

At ~ 10ms - 2s (B5) 
Consequently 
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At 



< At0 = T0, 



At 



(B6) 



'.max 



< 



Rc/C 



Rc/c{l - (3) 



where the relation < cos^^ < 1 has been apphed. Combining (Bl), (B5) and (B6), the 
range of Arg is hmited by 



For one typical value of T, namely r=100, then I - p = {T - VT^ - l)/r ~ 0.00005 is 
decided, therefore 



Under this approximation, Atq is allowed to take some representative values such as Are=0.001, 
0.01, 0.1, 1, 10 and so on. 
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